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Abstract 

Neutral pion photoproduction on the deuteron in the inelastic channel is 
studied in the threshold region from 142 to 160 MeV. The calculation is 
based on the use of the diagrammatic approach. Contributions from the pole 
diagrams, one-loop diagrams both with n-p and vr-A^ rescattering as well as 
two-loop diagrams with simultaneous inclusion of these rescattering mecha- 
nisms are taken into account. We have found a reasonable agreement with 
a very simple calculation by J.C. Bergstrom et al. [Phys. Rev. C 57 (1998) 
3203] . It is shown that the 'effective' electric dipole deuteron amplitude is 
negative in agreement with a prediction of ChPT. We conclude that a devia- 
tion of about 20% from a ChPT prediction for Eg^ discovered in an experiment 
on coherent vr*' photoproduction off the deuteron cannot be attributed to in- 
adequate estimates of inelastic channel contributions. We suggest that vr-A^ 
rescattering is responsible for the deviation of the free-nucleon pi threshold 
amplitude from the value measured in that experiment. 

PACS: 25.20.Lj; 25.45.De 
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I. INTRODUCTION 



Nowadays, the threshold s- and p- amplitudes of 7r° photoproduction off the deuteron 
are of considerable theoretical and experimental interest. This is motivated by two reasons. 
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Firstly, the deuteron is a natural 'source' of neutrons. If one wants to minimize uncertainties 
stemming from nuclear forces when extracting neutron threshold parameters one should 
use a deuteron target. At present, there is a significant disagreement for values of the 
threshold electric dipole amplitude E^^^ calculated in chiral perturbation theory (ChPT) 
and dispersion theory (DR), respectively, 2.13 and 1.19 [Q] (in units of 10~^/ ^1^,+ which are 
suppressed from here on). At the same time the predicted values for Eq^ agree reasonably 
with each other (—1.16 and —1.22, respectively) and with experimental values —1.32 ±0.08 
0] and —1.31 ± 0.08 0]. Resolving the problem mentioned in favour of the one or the 
other theory requires experimental measurements of which in the reactions 'yd — ^ ir'^d and 
ed —>■ e'TT^d seem to be most promising. However, even in these simplest reactions many 
theoretical problems in the extraction of neutron data are anticipated. 

On the other hand, the deuteron cannot be considered as a direct sum of the proton 
and neutron (impulse approximation). Photons may also interact with potential (virtual) 
mesons. This effect is known as 'Meson Exchange Currents' (MEC). It was found to be 
very important in the case of coherent neutral pion photoproduction on the deuteron 0. 
Accepting the ChPT predictions for free nucleons we could expect in the impulse approxi- 
mation a value of about +0.5 (= ^[E'o+ + -^0+"]) for the deuteron electric dipole amplitude 
Ed- In fact, however, the ChPT result is quite different being about —1.8 ± 0.2 and the 
difference is due to the meson exchange contributions. 

A very recent measurement of E^ performed at SAL in coherent 7r° photoproduction 
on the deuteron within 20 MeV of threshold has given a value of —1.45 ± 0.09 which is 
about 20% lower than the ChPT prediction. There is also a significant deviation from the 
ChPT prediction for the free-nucleon pi amplitude. A comparison of the differential cross 
sections calculated in a theoretical model with the data revealed a noticeable disagreement 
as well. It was claimed in Ref. that a possible reason for the disagreements might 
consist in inadequate treatment of the inelastic channel in Ref. 0. Since the ^i/(7, 7r°)np 
reaction cannot be resolved one has to estimate its contribution to the coherent channel 
making the use of theoretical predictions. Indeed, a model developed in Ref. seems to 
be oversimplified. For example, the authors employed a square-well deuteron wave function 
and used the effective range approximation to describe the final state n-p interaction. The 
deuteron d-wave was ignored. Explicit calculations of very important diagrams with vr-A^ 
rescattering were not carried out. Effectively these latter where included into the model by 
accepting the extracted value —1.45 for Ed which, as it has been noted, contains contributions 
from MEC or, in other words, from tt-N rescattering. Although a theoretical uncertainty of 
±25%, based on input parameters, was assigned in Ref. to the calculated cross sections 
it would be desirable to have more realistic calculations for the inelastic channel. 

The present paper is aimed at fixing the defects of the model . We are going to study 
whether these defects may be responsible for the deviations mentioned. In our previous 
paper P on the reaction 

■yd — s> TT^np (1-1) 

we restricted ourselves to the energy region 200 to 400 MeV. Now we extend the model 
to the threshold region. The extension includes: i) a more realistic treatment of the 7r° 
photoproduction amplitudes at threshold energies; ii) in Ref. P] we considered n-N rescat- 
tering mechanism but without taking into account the n-p rescattering process which may 
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follow Ti-N rescattering. In the present paper a corresponding diagram has been taken into 
account. We would like to emphasize that the study of the process ( p. .11) in the threshold 
region is hardly of theoretical interest in its own. Rather it is needed for experimentalists 
to estimate the relative contribution of the inelastic channel to the cross section of coherent 
pion photoproduction on the deuteron. Furthermore, these studies may be useful for extract- 
ing the elementary -Eo+° amplitude from exclusive inelastic data where the recoil neutron is 
detected. 



II. KINEMATICS 

Let us denote by = (tu, k), = (e^, -k), q = (e^, q), p„ = p„) and Pp = (Sp, p^) 
the 4-momenta of the initial photon and deuteron and the final pion, neutron and proton, 
respectively in the j-d cm. frame. A symbol we reserve for the lab photon energy 



= W^d^/fnd with W^d = uj + ed = uj + yuj'^ + m\ and being the deuteron mass. 
It is convenient to take as independent kinematical variables the photon energy and pion 
angle, Gtt, in the cm. frame and momentum p of one of the nucleons - say, the neutron - 
in the final n-p cm. frame (see Ref. 0). Once these variables have been specified the pion 
momentum can be easily found: 



|q| 



2W. 



'yd 



^[w',d - (w^np + lifW^.d - iW^P - i^Yl (2.1) 



where Wnp = 2-\/p^ + rri^ {m is the nucleon mass) and /j, is the 7r° mass. Then making a 
boost with the velocity — q/ {W^d — £n) we obtain the neutron momentum p„ and, therefore, 
totally restore the kinematics in the ^-d cm. frame. Note that in kinematic calculations 
a distinction is not made between the proton and neutron masses. It is, however, taken 
into account in the threshold 7r° photoproduction amplitude when parametrizing the s- and 
p- waves multipoles (see Sect. Eq. (p.2|)). 
The differential cross section is 

d'^a 1 2m^ed|qP 1 ^ i i, \i2 

"« \{mpmn\T\Xmd}\ , (2.2) 



dpd^l^ (27r)5 AujW^dWnpie^q ■ Pp - Sp^j?) 6 



d 



where rup, rUn, A, and rud are spin states of the proton, neutron, photon, and deuteron, 
respectively. 

Assuming that near threshold energies the nucleons are not detected we should integrate 
the l.h.s. of Eq. (^.21) over the momentum p to obtain 



da f d'^a ^ 



p^dpdQp, (2.3) 



dfl-j, J dpd^T, 

where the maximum value p'^'^^ is found from energy-momentum conservation and reads 



p-- = -^(H^^,-^)2_4m2. (2.4) 
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III. THEORETICAL MODEL FOR THE -fD tt^NP PROCESS 
The scattering amplitude of the reaction (|1 . ID is described as a sum of contributions from 



diagrams which are expected to be most important in the energy region under consideration. 
The pole diagrams a and b in Fig. |I] and one-loop diagram 0c correspond to the impulse 
approximation without and with n-p interaction in the final state, respectively. The necessity 
to consider the latter mechanism follows from the very strong nucleon-nucleon interaction 
at small energies. It might seem that diagrams to / with vr-A^ rescattering could be 
disregarded in the threshold region since the vr-A^ scattering lengths are about two orders 
smaller than those for n-p scattering. In fact, however, there is actually no suppression of 
the diagrams to / in comparison with diagram |l|c. Indeed, keeping in mind that the 
threshold electric dipole amplitudes i?o+ for charged channels are about 30 times larger in 
absolute numbers than those for the neutral channels one can expect the corresponding 
contributions from tc-N rescattering to be of the same order as the ones from final state n-p 
interaction. This expectation will be confirmed below by numerical calculations. A very 
big effect of tt-N rescattering in the case of coherent 7r° photoproduction on the deuteron 
in the threshold region was also found in Refs. Note that our treatment of vr-A^ 

contribution is similar to the one of 'three-body interactions to order g^' from Ref. in 
the case of coherent 7r° photoproduction off the deuteron. It was found in that paper that 
three-body contributions at order are less important. Based on this finding we do not 
take into account the corresponding diagrams for the reaction ( p..l[ ) though they of course 
exist. 

We begin with the pole diagrams. The matrix element corresponding to the diagram in 
Fig. |l]a reads 



(mpm„|T^''(k,q,Pp)|Amrf) = ^ (^Pp + 2 j (^n\T^n^^On{Kn,q.^n)\^mn), (3.1) 

where (p^ + k/2) is the deuteron wave function (DWF) and (m„|Ty„^^o„|Am^) is the 

amplitude of the elementary process 7^ —>■ n^n. The amplitude depends on photon (k^„) 
and pion (q^„) momenta taken in the cm. frame of the ir-n pair. These momenta can be 
obtained from the corresponding momenta in the 7-6? cm. frame through a boost with the 
velocity p^/ (W^d - £p)- 

In calculations we used DWF for non-relativistic versions of the Bonn OBE potential 
11,12| (the same models were used when solving the Lippmann-Schwinger equation for the 
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p scattering amplitude needed for calculations of diagrams |I|c and /, see below). Analytical 
parametrizations of the s- and d-wave amplitudes were taken from Ref. |13[. We would like 
to note just here that our results are practically independent of the choice of the potentials. 
In our previous paper P] we used the well-known Blomqvist-Laget parametrization 



of the pion photoproduction amplitudes. This model is rather simple and, nevertheless, 
fits reasonably to the experimental data. Unfortunately, it is obsolete in the sense that it 
predicts the threshold electric dipole amplitudes for the neutral channels in accordance with 
the classical low-energy theorem, Eq^ = —2.4 and -^0+° = 0.4, but in disagreement with the 
modern values (see Sect. 0). To be consistent with these latter, we use in the energy region 
uj < 150 MeV a parametrization for the neutral pion photoproduction amplitude given in 
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FIG. 1. Diagrammatic representation of the scattering amplitude. 

Ref. [^. The threshold amplitude in the cm. frame can be written in the following 
form []I| 

T^n^ttOn = - ^'^^^^ [ia- ■ e {Eo+ + k ■ q pi) + ia ■ k e ■ q p2 + q - (k x e) ps] , (3.2) 
m 

where W^n is the total energy of the system, and pi, p2 and ps are the p- waves multi- 
poles. The i?o+ amplitude was taken in the form of Eq. (25) of Ref. [0 with parameters from 
Table I of the same work. The parametrization takes into account the energy dependence 
of the £'0+ multipole as well as the so-called cusp effect. Just at threshold it reproduces the 
ChPT and DR values (see Sect. |). These parametrizations were used in the present paper 
to study uncertainties due to the choice of a theoretical model for pion photoproduction. As 
to the p-wave amplitudes, they are also taken for two models and listed in Table |I[ As was 
shown in Ref. , the parametrizations corresponding to ChPT and DR reproduce very well 
the available data on the energy dependence of the real part of i?o+ as well as the total and 
differential cross sections of the reaction 'jp —>■ ir^p up to ~ 180 MeV. 

One should emphasize that, although we are going to work in the energy region below 
160 MeV, we need parametrizations for the photoproduction amplitudes at higher energies 
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TABLE I. Threshold p- waves amphtudes (units are 10 ^/^^+). 







Pi 




P2 




TT^P 


P3 


TT^n 


ChPT 
DR 


10.3 
10.5 




7.4 
7.8 


-11.0 
-11.4 


-8.4 
-8.8 


11.7 

10.2 




11.1 ° 

9.5 " 



"The value was obtained from isospin relations between four physical channels. 

too. Such energies emerge because of two kinds of boosts. The first one occurs at the 
integration over p in Eq. ( |2.3| ). Below 160 MeV the value of p"^"-^ does not exceed 125 
MeV/c so that the corresponding boosts lead to only small energy shifts. But at integrations 
over loops when calculating n-p final state interaction and vr-A^ rescattering (see below), 
energy shifts can be much bigger. In addition, when considering vr-A^ rescattering the pion 
photoproduction amplitudes for the charged channels are also needed. Both these latter 
and the photoproduction amplitudes at photon cm. energies above 150 MeV were taken 
from the Blomqvist-Laget model fl^ . 

In full analogy with Eq. ( |3.1| ) we write the matrix element corresponding to diagram 
in the form 



[mpTTln 



l^^''(k,q,pJ|Amrf) = Jl^m'm^ (Pn + I) (mp|T^p^^Op(k^p,q^p)|Amp). (3.3) 

Tip \ ^ / 



The photon (k^^p) and pion (q^p) momenta in the cm. frame of the n-p pair can be obtained 
by boosting the corresponding momenta in the 'j-d cm. frame with the velocity Pn/C^jd — 

The one-loop diagram in Fig. |l]c with n-p rescattering in the final state is expected to 
be very important in the energy region under consideration. The smallness of the relative 
n-p momenta leads to strong ^Sq- and ^^i-wave interactions which significantly modify 
the n-p plane wave states appearing in the pole diagrams |l|a and h. The matrix element 
corresponding to diagram |I|c is 



(mpm„|T^''(k, q, p„)|Amd) = ~ m j 



X 

m,-,m 



{'^T^Y Pin - pint - 

H [{Pout^^P^n\Tnp\ - Pi„,mpmfi)(mpmfi|T^"(k,q,pJ|Amd) 



+ {Pout^^P^n\Tnp\Pin,mpmf,){mpmf,\T (k, q, pJ|Amd)J , (3.4) 

where p^^^ = (pp — p„)/2 and Pj„ = p^ + q/2 are the relative momenta of the n-p pair after 
and before scattering, respectively, and (p^^^, mpm„|T„p|p.^, mpm^^) being the half off-shell 
n-p scattering amplitude. We will not discuss here details of computations of the amplitude 
( p.4|) because they are given in Ref. [§. Note only that all partial waves with the total 
angular momentum J < 1 were retained in the n-p scattering amplitude. We cut off the 
integration in Eq. ( p.4| ) at |PsP"^ = 500 MeV/c since it is difficult to control the pion 
photoproduction amplitudes at higher momenta. Such a cutting off has no impact on our 
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final results because the integral is mainly saturated at smaller momenta. For example, a 
variation of |PsP"^ from 400 to 500 MeV changes the calculated differential cross sections 
by less than 2%. One further remark is concerned with the choice of particle energies when 



integrating over the momenta p^. In accordance with a prescription of Ref. [0 we will 
suppose a spectator particle (it is denoted in Fig. |l|c by a subscript s) to be on its mass 
shell. 

The amplitude for the diagram in Fig. reads 



mpmn\T (k, q,p„)|Amd) 



_ [ y I + 



(2vr) 




e* + e* 



2W, 



nn pf - Pin " ^0 



{nin I T^o„^^o„ (p;„ , p* ) I mi ) (mp I T^p^^Op (k^p, q- ) | ma ) 



■{mn\Ti' 



p^iTOn\fnni 



Pln,P*s)\^l){'mp\T^fn->ir-p{^np,(lip)\m2) 



(3.5) 



In Eq. 

et + e* 



151) the asterisk denotes kinematical variables in the ir^-n cm. frame. 



Pljf + /i^ + \ Pln^ + rn? is the total energy in this frame and 



Pn 



\Wl^ - (m + y.f\\Wl^ - (m - > 0. 



(3.6) 



is needed to transform the momentum p* to p^ 

rsi) 



A boost with the velocity — Pp/ (W^d — 

in the 7-6? cm. frame. As in the case of Eq. ( p.4|) we cut off the integration in Eq. 
at Ip*!™"^ = 500 MeV/c without any noticeable changes in the calculated differential cross 
sections. Spectator particles in diagrams and e are again supposed to be on their mass 
shells. 

The Ti-N scattering amplitude has the following form 



(m'|T;^_,^,(q',q)|m) = (2 



E 

JLm 1 



q)TM,q), (3.7) 



where q' and q are the final and initial momenta of (generally, off-shell) particles in the 
71-N cm. frame, Yj^^{q) are the spherical harmonics, Cj^_^j^^f^ are the Clebsch-Gordan 
coefficients. Partial wave amplitudes Tj'j^q', q) for all s- and p- waves channels were obtained 
by solving the Lippmann-Schwinger equation for a separable energy-dependent vr-iV potential 
built in Ref. [|1^]. Only for the psa-wave we preferred using a model from Ref. with 
explicit inclusion of a s-channel pole diagram with an intermediate A isobar. Contributions 
from the d-waves were found to be negligible. Two vr-iV scattering amplitudes in isospin 
space are numbered by a symbol / = | or |. Needed in Eq. (|3.5| ) physical amplitudes are 
linear combinations of these isospin amplitudes: 



T^0„^^0„ = I (T2 + 2T2 



(3.8) 
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In full analogy with Eq. ( p.5|) one has for diagram |l|e 



{mpmn\T ""(k, q, p )|Amd) 



J (2vr)3 i;;;, '"^"^^ \^ 



' 2 \ 2W^p pf-p*2_^Q 




{mp\T^Op^^Op{plp, p*) |mi) (m„|T^„^^o„(k^„, q-„) jma) 
-(mp|T^+„^^0p(p;p,p*)|mi)(m„|T^p^^+„(ks„,q-„)|m2) |, (3.9) 

where the meaning of all variables is easily understood from the corresponding ones in Eq. 

(PD- 

Finally, let us consider diagram |I|/ which includes simultaneously both tt-N and n-p 
rescattering. One has for the corresponding matrix element 



(mpm„|T^-^(k,q,p„)|Amd) = - m j 



X 



[{Vout^^P^n\Tnp\ - Pi„,mpmfi)(mpm^|T^'^(k,q,pJ|Amd) 



+ {^ouvT^v^n\Tnp\Vin^mpmn){mpmn\T (k, q, p J | Am^) J . (3.10) 

The amplitude ( p.lO| ) includes as subblocks the amplitudes ( p.5| ) and ( p.9| ) but with off-shell 
outgoing nucleons so that p*jv^ in Eq. ( p.6|) may be both positive and negative. One point 
which has to be emphasized is how to specify the energies in the loops. It is reasonable 
to suppose that the nucleon Ng in the left loop is on its mass shell since such a choice 
precisely corresponds to that used when evaluating diagrams ^c? and e. As to the right loop, 
a prescription for the nucleon energies used in the case of coherent 7r° photoproduction on the 
deuteron consisted in putting the recoil nucleon on its mass shell in the pion photoproduction 



vertex ||10[. This prescription (which is referred to from here on as Case 1) was motivated 
by the location of the leading singularities of elementary amplitudes (the A pole singularity) 
and internal particle propagators. It was found in Ref. |T^ that putting the recoil nucleon on 
its mass shell at the vr-A^ vertex led to a one and a half times larger rescattering amplitude. 
We think, however, that the prescription above can hardly be considered as a hard-and-fast 
rule. Even a glance on Fig. enables one to see that there are no reasons why such a choice 
of nucleon energies should be preferred. In fact, we have checked other realistic possibilities 
to specify the nucleon energies in the right loop. One of them was to suppose that each of 
the nucleons carried equal parts of the total energy of the final n-p pair, namely (e„ + ep)/2 



(Case 2). One more choice rejected in Ref. [|T0| was to put the recoil nucleon on its mass 
shell at the tt-N vertex (Case 3). 



Summations over polarizations of the particles in Eqs. ( p.l| )-( pTTOD as well as the three- 



dimensional integrations in Eqs. (|3.4| ), (|3.5|) and ( |3.9|) and six-dimensional one in Eq. ( |3.10|) 



have been carried out numerically making use of the methods described in Refs. p!7| , p!3 



Note, since Eq. (|2.3|) involves a three dimensional integration we were forced to carry out 
integrations with dimensions up to and including nine. Additional one-dimensional inte- 
grations were needed to solve the Lippmann-Schwinger equations for the A^-A^ and vr-A^ 
scattering amplitudes. The computing work was very hard but we were in a position with 
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reasonable computer resources to evaluate the integrals with a precision providing an accu- 
racy of calculated cross sections better than 10%. 



IV. RESULTS AND DISCUSSION 

Based on the approach described in the previous section, we calculated the differential 
cross section ( p^.3D and total cross section of the reaction ( |1.1| ) in the energy region from the 
threshold (E*'"' = 142.2 MeV) to 160 MeV. 




0.000 
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FIG. 2. Shown are contributions of different diagrams to the differential cross section ( |2.3| ) in 
the cm. frame. Contribution of the pole diagrams ||a and b is given in dotted Unes. Successive 
addition of the diagrams with n-p rescattering (|^c) , vr- rescattering (^d and e) and diagram [l[f 
gives the dashed, dash-dotted and sohd lines, respectively. The curves correspond to the ChPT set 
for the multipoles in Eq. ( |3.2| ) and Case 1 for nucleon energies in the right loop of diagram |^. 

In Fig. we present our predictions for the differential cross section with the ChPT set for 
the multipoles in Eq. and with the choice Case 1 for nucleon energies in the right loop 
of diagram |l|/'. One can see that the contributions of all the diagrams in Fig. |l] are of the same 
order so that none of the diagrams can be ignored in the calculation. The two pole diagrams 
give a forward-peaked angular distribution. This reflects the fact that the ChPT values of 
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the isospin averaged threshold electric dipole amplitude Eq'^^ = |[-Eo+ + -^0+°] — O-^ and 
the amplitude p'{^^ = + p"'""] = 8.9 are positive. Indeed, it is seen from Eq. ( p.2| ) 

that the amplitudes -£^0+^ and p^^^ are responsible for the forward-backward asymmetry in 
the differential cross section. An effect of n-p final state interaction (diagram |I|c) is strongly 
dependent on the photon energy. Just above threshold (e.g., at £"^=145 MeV) it leads to 
increasing the differential cross sections over the full angular region and this increase is 
attributed to the strong attractive n-p interaction in the singlet ^S'o-state. At higher photon 
energies the repulsive triplet '^S'l-interaction grows in importance. At E > 150 MeV the 
total effect of diagram |I]c consists in a noticeable lowering of the cross sections except in the 
backward angle region. 
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FIG. 3. Shown is the sensitivity of the differential cross section ( |2.3D in the cm. frame to 
different prescriptions for nucleon energies in the right loop of diagram [l|/: Case 1 (solid), Case 2 
(dashed), Case 3 (dotted). The curves correspond to the ChPT set for the multipoles in Eq. (|3.2D . 

As to diagrams and e with tt-N rescattering in the final state, they increase the 
differential cross section at all angles. At > 155 MeV these diagrams contribute mainly 
at backward angles. Finally, diagram |llf gives a noticeable reduction of the cross section at 
forward angles. It is seen in Fig. |^ that after inclusion of diagrams |l]d to / the differential 
cross section becomes to be backward peaked at E > 150 MeV. This means that due to 
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Tc-N rescattering the amplitude -^0+ effectively acquires a negative contribution altering its 
sign. When making such a statement we suppose (and with some justice, see below) that 
although rescattering can also modify the Pi^^ amplitude the latter remains to be positive. 
This effect was first discovered in the framework of the ChPT p in the study of coherent 
7r° photoproduction on the deuteron. 

With the ChPT set for the multipoles in Eq. ( |3.2| ) we have calculated and presented in 
Fig. ^ the differential cross sections for various choices of nucleon energies in the right loop of 
diagram One can see that the results at < 90° are practically independent of the choice 
of the energies. This is not the case for backward angles. Case 1 leads to minimum cross 
sections and Case 3 to maximum ones. Results with Case 2 lie between. Simultaneously, with 
increasing photon energy the deviation is declining. It is worth mentioning that analogous 
observations were made in Ref. in consideration of coherent vr" photoproduction in the 
threshold region. It was found in that paper that various approximations for the nucleon 
kinetic energy at integrations over loops gave noticeable changes in the calculated differential 
cross sections especially near threshold and backward angles. 
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FIG. 4. Shown is the sensitivity of the differential cross section ( |2.3D in the cm. frame to sets 
of the multipoles in Eq. ( |3.2D : DR (dashed), ChPT (solid). The curves correspond to Case 1 for 
nucleon energies. 
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In Fig. ^ we compare our predictions for the angular distribution of the differential cross 
section with the ChPT and DR sets for the multipoles. It is seen that the ChPT and DR 
sets produce rather close cross sections at backward angles for > 150 MeV though there 
exist noticeable deviations at forward angles. In the energy region under consideration the 
DR sets produce only half of the ChPT cross sections at forward angles. This reflects the 
significant difference between the ChPT and DR predictions for the threshold electric dipole 
amplitude -Eo-f" (^^e Sect. |) and the close consistency of the two for the p{^^ amplitude. 

The previous consideration shows that t^-N rescattering (diagrams ^d-f) markedly 
changes the results calculated without this effect (diagrams |I|a-c). It is instructive to study 
the question: what modification of the low energy parameters in the model without t^-N 
rescattering is required in order to reproduce the predictions of the full model. We will be 
interested in the amplitudes E^y^ and p\^^ only. Let us suppose that due to vr-A^ rescattering 
these amplitudes acquire additional contributions ^e'^ and ^p\^\ Numerical estimates 
for them are obtained by switching off diagrams |l|d-/and adjusting these parameters to re- 
produce the differential cross sections at 9,^ = 0° and 180° for the full model. Our results are 
given in Table O. One comment should be made here. In a rigorous consideration the values 
^Eq^ and ^pi^ must be independent of the parametrization for the low energy amplitudes 
in Eq. ( ^.2| ), i.e. the former should be the same for ChPT and DR amplitudes. Indeed, the 
difference between ChPT and DR manifests itself in the vr'' production amplitudes but the 
main contribution from tt-N rescattering stems from charged pion exchanges. The slight dif- 
ference for ChPT and DR results in Table |I| is because of a rather rough numerical method 
of extracting ^Eq^ and ^Pi'\ Nevertheless, we beheve that such a method properly repro- 
duces the tendency of variation of E^^ and pj^"* due to vr-A^ rescattering. The dependence 
of ^Eq^ and Ap^"'"'' on the choice of nucleon energies as seen in Table |l| is expected since 
the Cases 1 to 3 effectively lead to different pion propagators in Eq. ( |3.10| ). 

Averaging over all values in Table || we obtain /S.Eq^ ^ —1.1. One can see that A^o? 
is negative in accordance with the ChPT result from Ref. ^ being, however, nearly two 
times smaller in the absolute value. The p^^"* amplitude acquires a positive contribution an 
averaged value for which is 4.1 ± 0.6. This result being put together with the free-nucleon 
value of about 9 gives pi = 13.1 ± 0.6 and reasonably reproduces an experimental result 
of Ref. [0, pI^^ = 12.88 ± 0.28. In other words, our estimates confirm an assumption from 
that paper that a mechanism responsible for renormalizing the pi amplitude is the two-body 
process corresponding to vr-A^ rescattering. 



TABLE II. Effective modifications of the threshold Eq_^' and p\ amplitudes due to vr-A^ 
rescattering. Units are conventional. 





Case 1 




Case 2 




Case 3 












ApS+) 




ApS+) 


ChPT 


-1.2 


+3.7 


-1.1 


+4.2 


-1.2 


+4.7 


DR 


-1.0 


+3.5 


-1.1 


+3.8 


-1.1 


+4.5 



To make it much easier using our results in practice, we have parametrized the differential 
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cross sections at four photon energies 145, 150, 155, and 160 MeV with the help of the 
following formula 



^'^ Ai + A2Cose^, (4.1) 



with parameters given in Table [m| . The values for Ai and A2 have been obtained by a fit to 
the average of the minimum and the maximum cross sections (see above) and the errors in Ai 
account for variations of the cross sections both due to various choices of sets for multipoles 
and nucleon energies in the right loop of diagram |l]/. The corresponding cross sections at 
arbitrary energies up to 160 MeV can be obtained by making use of a suitable interpolation 
procedure. We think that such a simple parametrization provides reasonable accuracy to be 
used in estimates of inelastic channel contributions to coherent 7r° production. 



TABLE III. Coefficients Ai and A2 in Eq. (4J.) at four selected photon energies. Units are 
/ib/sr. 





145 MeV 


150 MeV 


155 MeV 


160 MeV 


Ai 


0.0043 ± 0.0014 


0.023 ± 0.005 


0.057 ± 0.007 


0.095 ± 0.012 


A2 


-0.0007 


-0.007 


-0.026 


-0.041 



2.5 




FIG. 5. Total cross section of the reaction '^H{'y,7r^)np in the threshold region as a function 
of photon lab energy. Meaning of curves as in Fig. ^. 

Performing in Eq. ( |2.3| ) an integration over the pion solid angle we can calculate the 
total cross section atot of the reaction ( |1.1| ) (note that atot = 'inAi where cxtot is in units of 
fih). In Fig. ^ we show contributions of different diagrams to o"tot with the ChPT set for 
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the multipoles in Eq. (|3.2| ). The two pole diagrams produce a total cross section which is 
rapidly increasing with photon energy. After inclusion of the n-p final state interaction the 
energy distribution of atot becomes to be more flat. Just above threshold this interaction in 
the ^S'o-wave leads to an increase of atot but at > 155 MeV the cross section diminishes 
by a factor of about 2. Diagrams |I]g? and e with ti-N rescattering give noticeable positive 
contribution to the total cross section. Of all the diagrams is the least important one, 
decreasing atot by about 10%. 




FIG. 6. Left panel: shown is the sensitivity of the total cross section to different prescriptions 
for nucleon energies in the right loop of diagram ||/. Meaning of curves as in Fig. ^. Right panel: 
sensitivity of the total cross section to sets of the multipoles in Eq. ( p.2[ ) . Meaning of curves as in 
Fig. |. 

The total cross sections for various choices of nucleon energies in the right loop of diagram 
|l]/ are presented in the left panel of Fig. p. It is seen that they are practically independent 
of a prescription for the energy. This reflects the smallness of the diagram contribution to 
the total cross section. Our predictions for atot with the ChPT and DR sets of the multipoles 
are also shown in Fig. There is a deviation of about 15% in the full energy region. 

In Fig. we compare our results with the ones of the model for the differential and 
total cross sections. One can see that the predictions of the two calculations are quite 
consistent. There is some disagreement at > 155 MeV and forward angles for the 
differential cross section. It should be emphasized, however, that this disagreement cannot 
make any noticeable changes in the cross sections of coherent 7r° photoproduction measured 
in Ref. because of smallness of the inelastic channel cross sections in comparison with 
the coherent channel ones. At the same time, at backward angles the inelastic channel is 
growing in importance but there the deviation is not very significant. One can see that the 
maximum deviation between the present calculation and the one of Ref. does not exceed 
20% at 160 MeV and backward angles, being even smaller when the energy is decreasing. 
(Note here that the authors of that paper ascribed to their calculated cross sections of the 
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FIG. 7. Left panel: shown are errors bars (filled areas) in the differential cross sections 
stemming from various prescriptions for nucleon energies and different sets of multipoles at 150, 
155, and 160 MeV (from bottom to top) calculated with Eq. ( ^.Ij ) and coefficients Ai and A2 from 
Table ID. Curves are results of Ref. Right panel: the same as for the left panel but for the 



total cross section. 



inelastic channel uncertainties of about ±25% which were taken into account in their data 
analysis). In the case of the total cross section we observe good agreement with Ref. 0. 



An analysis performed by Bergstrom [|18[ has shown that the pi amplitude is insensitive to 
the variation in the inelastic cross section from the present work and that from Ref. [Q. 
Likewise, the extrapolation of the Ed amplitude to threshold changes very little, remaining 
within the error in Eq. (12) of Ref. ||^. 

In summary, we have performed an investigation of the reaction ^i^(7, Tc^)np in the 
threshold region and found a good agreement with the only previous calculation for the 
total cross section. A reasonable agreement also exists for the differential cross sections. 
Our predictions coincide within of 20%. If one takes into account that the authors of 
that work ascribed to their results an uncertainty of about ±25% we conclude that the 
present calculation cannot lead to any changes in the total and differential cross section of 
coherent 7r° photoproduction on the deuteron measured in Ref. |Q. This means that the 
deviation of about 20% from theoretical predictions found in that paper for the threshold 
amplitude still remains to be resolved. Our calculations have shown that the difference of 
about 4 ■ 10~^//i^+ between the free-nucleon value for the pi threshold amplitude and the 
experimental value can be attributed to two-body effects due to t^-N rescattering. 
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